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We perform a relativistic chiral effective-field theory calculation of pion electroproduction off 
the nucleon (e~ N — > e~ N n) in the A(1232)-resonance region. After fixing the three low-energy 
constants, corresponding to the magnetic (Ml), electric (E2), and Coulomb (C2) yNA couplings, 
our calculation provides a prediction for the momentum-transfer and pion-mass dependence of the 
'yNA form factors. The prediction for the pion-mass dependence resolves the discrepancy between 
the recent lattice QCD results and the experimental value for the "C2/M1 ratio" at low Q 2 . 
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The A(1232)-resonance, the first excited state of the 
nucleon, dominates many nuclear phenomena at energies 
between the one- and two-pion production thresholds. 
The electromagnetic excitation of the A-resonance, the 
-fNA transition, has recently received a lot of attention. 
At low momentum-transfer (Q 2 ) it highlights the role of 
the pion cloud 00HSSlSl3> wnereas at larger j 
probes the onset of the perturbative QCD regime 

The jNA transition is predominantly of the magnetic 
dipole (Ml) type which, in a simple quark-model picture, 
is described by a spin flip of a quark in the s-wave state. 
Any ci-wave admixture in the nucleon or the A wave- 
functions allows for the electric- (E2) and Coulomb- (C2) 
quadrupole transitions. Therefore by measuring these 
one is able to assess the presence of the d-wave compo- 
nents and hence quantify to which extent the nucleon or 
the A wave-function deviates from the spherical shape 
( "hadron deformation" ) [ljj ■ 

The 7./VA transition has been accurately measured in 
the pion photo- and electro-production reactions Q, 0, 
0, E} • The E2 and C2 are found to be relatively small, 
the ratios Rem = E2/M1 and R S m = C2/M1 are at 
the level of a few percent. On the theoretical s ide, the 
most recent state-of-the-art lattice QCD study ob- 
tained a puzzling result: the computed ratio Rsm at low 
momentum-transfer appears to be significantly different 
from the observed value, ft is important to note that 
the lattice calculations were done at larger pion masses, 
while the result compared with experiment was obtained 
by a linear extrapolation to the physical pion mass. 

In this Letter we present a first chiral effective-field 
theory (%EFT) calculation of pion photo- and electro- 
production on the nucleon in the A-resonance region. 
Besides finding a good agreement of our calculation with 
observables, we are able to study the chiral behavior (m^- 
dependence) of the 7iVA transition. Our results show 
that there is no apparent discrepancy between the lat- 
tice data and the experimental result for Rsm- 

Our starting point is the relativistic chiral Lagrangian 
of pion and nucleon fiel ds 1121 supplemented with the rel- 
ativistic A-isobar fields |L3j . We organize the Lagrangian 



D-^ , such that superscript i stands for the power of elec- 
tromagnetic coupling e plus the number of derivatives 
of pion and photon fields. Writing here only the rele- 
vant terms involving the spin-3/2 isospin-3/2 field ij) 1 * of 
the A-isobar we have (with antisymmetric products of 
7-matrices: Y v = §[7^,71, = ie^^lplb)- 
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where M ~ 0.939 and Ma — 1.232 GeV are, respectively, 
the nucleon and A-isobar masses, TV and ir a (a — 1,2, 3) 
stand for the nucleon and pion fields, is the covariant 
derivative ensuring the electromagnetic gauge-invariance, 
and P" are the electromagnetic field strength and 
its dual, T a are the isospin 1/2 to 3/2 transition matrices, 
/tt ~ 92.4 MeV is the pion decay constant. C A contains 
the Rarita-Schwinger Lagrangian 0] of a free spin-3/2 
field formulated such that the number of spin degrees 
of freedom is constrained to the physical number. The 
couplings in Eq. are consistent with these constraints 
because of a spin-3/2 gauge symmetry [li). 

We next turn to the power-counting for the pion 
electroproduction amplitude using the "^-expansion" 
scheme 0|. In this scheme the excitation energy of the 
A-resonance: A = Ma — M ~ 0.3 GeV is treated as 
a light scale, so that for A ~ 1 GeV representing the 
heavy scales in the theory, we can use a small parameter 
S = A/ A. The other typical light scale of the theory, the 
pion mass, is counted as two powers of the small param- 
eter: TOtt/A ~ S 2 . The latter rule is the main distinction 
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sider the general form of the jNA vertex: 
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u a (p') r; NA u( P ) -- V2M [(Ma + M)2 + Q2] 

x iZ Q (p') {3a/(Q 2 )£^ kA pUa 

+ 9E(Q 2 )(q a p^-q-p'g a n^o (2) 

+ 3c(Q 2 )(g Q 9 M -9 2 .9 QM )^ 5 }"(p), 



FIG. 1: Diagrams for the eA^ — > eirN reaction at NLO in the 
5-expansion, considered in this work. Double lines represent 
the A propagators. The crossed nucleon-exchange graph is 
not shown in (b), but is included in the calculation. 



of this scheme from the previous power countings |lfill7l | 
which count A and at the same order. This difference 
plays a crucial role in separating the low-energy and res- 
onance regimes, as well as in approaching the chiral limit 
where m, v vanishes while A remains finite. Because of the 
distinction of m n and A the counting of a given diagram 
depends on whether the characteristic momentum p is in 
the low-energy region (p ~ m n ) or in the resonance region 
(p ~ A). In the resonance region, one distinguishes the 
one- A- reducible (OAR) graphs 0], see e.g., graph (a) in 
Fig. n Such graphs contain A propagators which go as 
l/(p — A) and hence forp ~ A they are large and all need 
to be included. Their resummation amounts to dressing 
the A propagators so that they behave as l/(p — A — S) . 
The self-energy S begins at order p 3 and thus a dressed 
OAR propagator counts as l/<5 3 . 

The pion electroproduction amplitude to next-to- 
leading order (NLO) in the <5 expansion, in the reso- 
nance region, is thus given by graphs in Fig. [J a) and 
(b), where the shaded blobs in graph (a) include correc- 
tions depicted in Fig.^c-f). The hadronic part of graph 
(a) begins at O(S ) which here is the leading order. The 
Born graphs Fig. ^b) contribute at 0(5). We note that 
at NLO there are also vertex corrections of the type (e) 
and (f) with nucleon propagators in the loop replaced by 
the A-propagators. However, adopting the on-mass shell 
renormalizations and Q 2 <C AA, these graphs start to 
contribute at next-next-to-leading order (NNLO). 

We have not shown the 7 AT A- vertex correction graph 
where the photon couples into the ttNN vertex, because 
at this order the effect of this graph can fully be absorbed 
in the graphs Fig.QJe) and (f) by a field redifmition relat- 
ing the pseudovector and pseudoscalar 7rAW couplings. 
Having done that, we compute graphs Fig. ^e) and (f) 
using the pseudoscalar coupling. 

The self-energy correction, Fig. ^c), was computed 
previously I n that calculation, the experimental 

value for the A-resonance width fixes h A — 2.85. To 
present the results for the vertex corrections we first con- 



where u a is the A vector-spinor, u is the nucleon spinor, 
q = p' — p is the photon 4- momentum, Q 2 = —q 2 , and 
9m i gE, and gc are the form factors which at Q 2 — are 
equal to the physical values of corresponding parameters 
from Lagrangian Q. These form factors relate to the 
conventional magnetic {G* M ), electric (G %) a nd Coulomb 



(G* c ) form factors of Jones and Scadron [ljj as follows: 
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where Q± = ^{M A ±M) 2 + Q 2 , Q 2 = Q 2 /Ml, /3 7 = 



1(1 



Q 2 ), with r = M/M A . The ratios E2/M1 



and C2/M1 at the resonance position can be expressed 
in terms of these form factors as: 
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The one-loop corrections to the 7iVA form factors are 
given by the graphs in Fig. ^e) and (f). For example, 
the (MS-subtracted) result for the graph (e) in Fig. 
can be cast in the form: 
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1 1-2/ 

-C NA J dyy(2y - 1) J dx 

o o 
x [xr + (1 - x - y) (1 + r)] M~ 2 , 

where M 2 = (x - (3) 2 - A 2 + 2(5^xy + Q 2 y(l - y) - is, 
fi = mjM A , /3 = 1(1 - r 2 + M 2 ), X 2 = /3 2 - M 2 , 
Cna = 4,g A h A Q 2 + /[3{l + r){87rU) 2 ], g A ~ 1.26. Anal- 
ogous expressions are obtained for the graph Fig. df). 
Alternatively, we have computed these graphs by using 
the sideways dispersion relations (see, e.g., |20j), and ob- 
tained identical results. 



cos 



cos 



FIG. 2: (Color online) xEFT NLO results for the 6^ depen- 
dence of the 7*p — > 7T p cross sections at yfs = 1.232 GeV and 
Q 2 = 0.127 GeV 2 . The theoretical error bands are described 
in the text. Data points are from BATES experiments pll2l|. 



The vector- meson diagram, Fig. dd), contributes to 
NLO for Q 2 ~ AA.We include it effectively by giving 
the <7M-term a dipole Q 2 -dependence (in analogy to how 
it is usually done for the nucleon isovector form factor): 
9m -> 9m (1 + Q 2 /0.71 GcV 2 )~ 2 . Analogous effect for the 
gE and gc couplings begins at NNLO and is not included 
in the present calculation. 

We now present the electroproduction observables cor- 
responding to the NLO amplitude of Fig. ^ Denoting 
the invariant mass of the final ttN system by s, we re- 
strict ourselves to the resonance kinematics: s = M\. 
The 7*iV — > irN cross section for unpolarized nucleons 
are expressed in terms of 5 response functions as : 



da 



dan 



dar 



-\/2e(l + e) cos$ 



dQ-K 
daLT 



s cos 2$ ■ 



da' 



TV 



dfln 



(6) 



hy/2e(l 



. , da'rrp 



where O w and $ are the pion polar and azimuthal cm. 
angles, respectively, and h denotes the electron helicity. 

In Fig. we show our xEFT results for the different 
cross sections entering Eq. (0. The only free parame- 
ters in this calculation are the low-energy constants from 
Eq. Q , which were chosen to yield the best description of 
the data as g M = 2.88, g E = -1.04, g c = -2.36. Within 
xEFT, we can estimate the theoretical uncertainty of the 
NLO result due to higher-order effects. The NNLO cor- 
rections to the amplitudes are expected to be of order of 
5 2 , to^/A, or Q 2 /A 2 . Therefore, the theoretical uncer- 
tainty R err of an observable R, which involves a product 
of two amplitudes, is estimated as (taking here A = M): 
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FIG. 3: (Color online) Q 2 dependence of the NLO results for 
Rem (upper panel) and Rsm (lower panel). The blue circles 
are data points from MAMI for Rem Q , and R S m [2 l l23| | . 
The green squares are data points from BATES 0. 



where R av is an average value of R. In Fig.|21the average 
is taken over the range of 9 T . One sees that the NLO 
^EFT calculation, within its accuracy, is consistent with 
the experimental data for these observables. 

In Fig.[3]we show the Q 2 dependence of the ratios Rem 
and Rsm- Having fixed the low energy constants gu, gE 
and gc, the Q 2 dependence follows as a prediction. The 
theoretical uncertainty here (shown by the error bands) 
is estimated according to Eq. J7J with the average R av 
taken over the range of Q 2 from to 0.2 GeV 2 . From the 
figure one sees that the NLO calculations are consistent 
with the experimental data for both of the ratios. 

In Fig. 0] we show the m T dependence of the 7^A 
transition ratios, with the theoretical uncertainty esti- 
mated according to Eq. Q where R av is taken over the 
range of to 2 from to 0.15 GeV 2 . The study of the m w 
dependence is crucial to connect to the lattice QCD re- 
sults, which at present can only be obtained for larger 
pion masses. The recent state-of-the-art lattice calcula- 
tions of these ratios use a linear, in the quark mass 
(m q oc to 2 ), extrapolation to the physical point, thus as- 
suming that the non-analytic m 9 -dependencies are neg- 
ligible. The thus obtained value for Rsm at the physical 
run value displays a large discrepancy with the experi- 
mental result, as seen in Fig.^J However, our calculation 
demonstrates that the non-analytic dependencies are not 
negligible. While at larger values of m„, where the A is 
stable, the ratios display a smooth to^ dependence, at 
TOtt = A there is an inflection point, and for < A 
the non-analytic effects are crucial, as was also observed 
for the A-resonance magnetic moment 0] • The 
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work presented here will, 
efforts. 



hopefully, complement these 
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FIG. 4: (Color online) m w dependence of the NLO results 
at Q 2 =0.1 GeV 2 for Rem (upper panel) and Rsm (lower 
panel). The blue circle is a data point from MAMI |22^. 
the green squares are data points from BATES 0. The 
solid black diamonds are lattice calculations |Tl|. whereas the 
dashed lines and open diamonds represent their extrapolation 
assuming linear dependence in m 2 . 



dependence obtained in xEFT clearly shows that the lat- 
tice results for Rsm may in fact be consistent with ex- 
periment. 

In conclusion, we have performed a manifestly gauge- 
and Lorentz-invariant %EFT calculation of the eN —> 
eNn reaction in the A(1232) resonance region. To NLO 
in the ^-expansion, the only free parameters entering the 
calculation are the jNA couplings gM, 9e, gc charac- 
terizing the Ml, E2, and C2 transitions. Our results 
agree well with recent high-precision data from MAMI 
and MIT-Bates at low Q 2 . The xEFT framework plays 
a dual role, in that it allows for an extraction of reso- 
nance parameters from observables and predicts their m w 
dependence. In this way it may provide a crucial connec- 
tion of present lattice QCD results obtained at unphysical 
values of m„ to the experiment. We have found that the 
opening of the A — > 7rA^ decay channel at m n — Ma — M 
induces a pronounced non- analytic behavior of the Rem 
and Rsm ratios. While the linearly-extrapolated lattice 
QCD results for Rsm are in disagreement with experi- 
mental data, the ;\EFT prediction of the non-analytic de- 
pendencies has allowed us to reconcile these results with 
experiment. As the next-generation lattice calculations 
of these quantities are on the way [25| , the %EFT frame- 
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